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Description 



[A METHOD FOR CONSIDERING 
HIERARCHICAL PREEMPTIVE DEMAND 
PRIORITIES IN A SUPPLY CHAIN 
OPTIMIZATION MODEL] 

Cross Reference to Related Applications 

[0001] The present application is related to pending U.S. Patent 
Application 10/ , filed concurrently herewith to Den- 
ton et al., entitled " A METHOD FOR SUPPLY CHAIN COM- 
PRESSION" having (IBM) Docket No. BUR920030197US1; 
U.S. Patent Application 10/ , filed concurrently here- 
with to Denton et al., entitled "METHOD FOR PURCHASE 
ORDER RESCHEDULING IN A LINEAR PROGRAM" having 
(IBM) Docket No. BUR92004009US1; U.S. Patent Applica- 
tion 10/ , filed concurrently herewith to Denton et 

al., entitled " A METHOD FOR SUPPLY CHAIN DECOMPOSI- 
TION" having (IBM) Docket No. BUR920040007US1; U.S. 
Patent Application 10/ , filed concurrently herewith 



to Denton et al., entitled "A METHOD FOR OPTIMIZING 
FOUNDRY CAPACITY" having (IBM) Docket No. 

BUR920030195US1; U.S. Patent Application 10/ , 

filed concurrently herewith to Denton et al., entitled 
"METHOD FOR FAIR SHARING LIMITED RESOURCES BE- 
TWEEN MULTIPLE CUSTOMERS" having (IBM) Docket No. 

BUR920040010US1; U.S. Patent Application 10/ , 

filed concurrently herewith to Denton et al., entitled " 
METHOD FOR SIMULTANEOUSLY CONSIDERING CUSTOMER 
COMMIT DATES AND CUSTOMER REQUEST DATES" having 
(IBM) Docket No. BUR920040008US1; and U.S. Patent Ap- 
plication 10/ , filed concurrently herewith to Orzell 

et al., entitled "METHOD FOR IDENTIFYING PRODUCT AS- 
SETS IN A SUPPLY CHAIN USED TO SATISFY MULTIPLE CUS- 
TOMER DEMANDS" having Docket No. BUR820030346US1. 
The foregoing applications are assigned to the present 
assignee, and are all incorporated herein by reference. 
Background of Invention 

[0002] piQi^ of Invention 

[0003] The present invention relates to computer implementable 
decision support systems for determining solutions to 
supply chain optimization problems in which customer 



demands have associated preemptive priorities to denote 
tlie relative importance of satisfying tlie demand. General 
methodologies within this field of study include advanced 
planning systems, optimization and heuristic based algo- 
rithms, constraint based programming, and simulation. 
[0004] Description of the Related Art 

[0005] A fundamental problem faced in all manufacturing indus- 
tries is the allocation of material and capacity assets to 
meet end customer demand. Production lead times neces- 
sitate the advance planning of production starts, inter- 
plant shipments, and material substitutions throughout 
the supply chain so that these decisions are coordinated 
with the end customers" demand for any of a wide range 
of finished products (typically on the order of thousands 
in semiconductor manufacturing). Such advance planning 
depends upon the availability of finite resources which in- 
clude: finished goods inventory, work in process inventory 
(WIP) at various stages of the manufacturing system, and 
work-center capacity. Often, there are alternative possi- 
bilities for satisfying the demand. Products may be built at 
alternative locations and within a location there may be 
choices as to which materials or capacity to use to build 
the product. The product may be built directly or acquired 



through material substitution or purchase. When limited 
resources prevent the satisfaction of all demands, deci- 
sions need to be made as to which demand to satisfy and 
how to satisfy it. This resource allocation problem is often 
addressed through linear programming. 
[0006] The be low- referenced U.S. Patents disclose embodiments 
that were satisfactory for the purposes for which they 
were intended. The disclosures of both the below- 
referenced prior U.S. Patents, in their entireties, are 
hereby expressly incorporated by reference into the 
present invention for purposes including, but not limited 
to, indicating the background of the present invention and 
illustrating the state of the art: U.S. Patent 5,971,585, 
"Best can do matching of assets with demand in micro- 
electronics manufacturing," October 26, 1999; U.S. Patent 
5,943,484, "Advanced material requirements planning in 
microelectronics manufacturing," August 24, 1999; and 
Nemhauser, G.L and Wolsey, LA., 1999, Wiley, Integer and 
Combinatorial Optimization. Also incorporated for reference 
into the present invention for purposes including, but not 
limited to, indicating the background of the present in- 
vention and illustrating the state of the art: (Leachman, R. 
C, R. F. Benson, C. Liu, D.J. Raar. 1996. IMPReSS: An au- 



tomated production-planning and delivery-quotation sys- 
tem at Harris Corporation-semiconductor sector, interfaces 
26(1) 637.) 
Summary of Invention 

[0007] The invention provides a method and system for deter- 
mining a production plan that allocates resources to dif- 
ferent demands having priorities by iteratively solving 
mathematical linear programs. Each mathematical linear 
program optimizes according to one of a plurality of sets 
of priorities and each iterative solution is consistent with 
the previous set of priorities. These priorities are hierar- 
chical and comprise two or more levels of hierarchy. 

[0008] Backorder costs penalties are determined independently 
for each set of priorities and comprise a full spectrum 
range within each set of priorities. The mathematical lin- 
ear programs solved in each iteration use the solution to 
the previous mathematical linear program as a starting 
point. The invention also adds constraints to the mathe- 
matical linear programs at each iteration to ensure that 
solutions to subsequent iterations are consistent with 
previous priorities. This method uses a different mathe- 
matical linear program for each iteration. The invention 
solves the mathematical linear programs for higher priori- 



ties before solving for lower priorities. 

[0009] More specifically, the invention provides a method of allo- 
cating resources to a hierarchy of demand priorities in a 
linear programming production planning system. In par- 
ticular, the invention aggregates the demand priorities 
into different priority groups and allocates the resources 
to the highest priority group of demand priorities using a 
first linear programming model. Next, the invention allo- 
cates remaining resources to the next highest priority 
group of demand priorities using a second linear pro- 
gramming model. The second linear programming model 
uses results from the first linear programming model. The 
invention continues this process by iteratively repeating 
the process of allocating remaining resources to the re- 
maining groups of demand priorities, in order of priority. 

[0010] When repeating the process of allocating remaining re- 
sources, the method uses a different linear programming 
model for each iteration and each different linear pro- 
gramming model uses results of the previous linear pro- 
gramming model. During the allocating processes, each 
linear programming model fixes variables for all demand 
priorities except for demand priorities of the priority 
group to which the resources are currently being alio- 



cated. Also, during the allocating processes, each linear 
programming model allocates the range of backorder 
costs within the priority group to which the resources are 
currently being allocated. The invention can also divide 
the priority groups into different sub-priority tiers. These 
sub-priority tiers can be processed simultaneously or 
separately. 

[0011] The present invention is a hybrid method which handles 
hierarchical demand priorities via demand class aggrega- 
tion and subsequent demand tier disaggregation.The in- 
ventive approach is based on aggregating demand priori- 
ties into groups, in which each group corresponds to an 
LP run. The invention generates backorder cost penalties 
associated with demand priorities within a group to model 
the consideration of multiple demand classes within a de- 
mand class group in the LP model. Further, the invention 
method iteratively modifies and solves the respective LPs 
in memory and leverages the results from one solution to 
the next. 

[0012] The inventive method results in substantial computational 
savings over conventional systems while at the same time 
respecting demand classes throughout the supply chain. 
The invention solves relaxed (e.g., modified) versions of 



the LP model at each iteration that allow for flexibility in 
realigning resources (e.g. work-center capacity) to ac- 
commodate lower priority groups without sacrificing 
higher priority groups. Further, the invention can be ap- 
plied to hierarchical demand priorities by selectively dis- 
aggregating demand tiers associated with demand 
classes. 

Brief Description of Drawings 

[0013] Figure 1: Overview of the structure of a linear program- 
ming application. 

[0014] Figure 2: Illustration of backorder coefficients within a 
priority group. 

[0015] Figure 3: Illustration of demand priority aggregation. 

[0016] Figure 4: Illustration of hierarchical demand priority based 

on demand class and demand tier. 
[0017] Figure 5: Illustration of the steps of the method. 

[0018] Figure 6: Flowchart illustrating one embodiment of the in- 
vention. 
Detailed Description 

[0019] Many factors are used to determine the relative impor- 
tance of the variety of customer demands that are ser- 
viced by a manufacturing firm (e.g. gross margin, strate- 



gic importance, forecast vs. committed orders, demand 
mix considerations etc.). Typically, the implementation of 
supply-chain planning methods requires that these be 
translated into a "demand priority" (e.g. 1, 2, 3,...) so that 
the many customer demands (measured in tens of thou- 
sands for division central runs) can be rank ordered ac- 
cording to importance. It is necessary to model this rela- 
tive importance of demands when trading off the alloca- 
tion of limited resources. For instance, photolithography 
tools are typically a production bottleneck in semiconduc- 
tor manufacturing. When insufficient capacity exists to 
schedule production for all demands it is expected that 
work-in-process (WIP) inventory associated with higher 
priority demands will be processed before competing 
lower priority WIP. 
[0020] It is important to break out the demand priority into a hi- 
erarchical set of attributes. For instance, it is common for 
a manufacturing firm to have two attributes a) demand 
class and b) demand tier which are used to determine the 
ultimate demand priority. Figure 4 illustrates this hierar- 
chical structure with an example in which there are four 
demand classes and each of the demand classes has three 
associated tiers. Thus, demand class is the most impor- 



tant factor and, within a particular demand class, impor- 
tance is dictated by demand tier. Tlie significance of tlie 
liierarcliical structure is to denote tliat demand class is 
the most important factor dictating priority while demand 
tier is secondary. 

[0021] To contrast the present invention, a conventional produc- 
tion planning linear program "LP" is shown below (such as 
that described in U.S. Patent 5,971,585, which is incorpo- 
rated herein by reference). This LP makes decisions in- 
cluding: production starts, material substitutions, and 
shipments planned to customers, between manufacturing 
and distribution locations, and from vendor suppliers. A 
LP is composed of an objective function that defines a 
measure of the quality of a given solution, and a set of 
linear constraints. The types of equations used in produc- 
tion planning models are well know to those practiced in 
the art and include:(l) Material Balance Constraints, which 
ensure conservation of material flow through the network 
of stocking points comprising the supply chain. 

[0022] (2) Capacity Constraints, which ensure that the capacity 
available for manufacturing activities is not exceeded. 

[0023] (3) Backorder Conservation Constraints, which balance the 
quantity of a given part backordered in a given planning 



period with tlie quantity bacl<ordered in the previous 
planning period and the net of new demand and new 
shipments. 

[0024] (4) Sourcing Constraints, which define target ranges 

(minimum and maximum) of shipments that should be 
made from a particular manufacturing or vendor location 
in the supply chain. 

[0025] A conventional LP formulation is provided below in the 
form familiar to those practiced in the art; i.e., definition 
of subscripts, definition of objective function coefficients, 
definition of constants, definition of decision variables, LP 
formulation or equations. 

[0026] Definition of Subscripts 

[0027] j- time period 

[0028] - material (part number) 

[0029] a - plant location within the enterprise 

[0030] „ _ material being substituted 

[0031] z _ group (which represents a family or collection of part 
numbers) 

[0032] g _ process (a method of purchasing or manufacturing a 
material at a plant) 



[0033] y_ receiving plant location 

[0034] _ demand center (i.e., customer location) (Note: the set 
of customer locations is mutually exclusive from the set of 
plant locations) 

[0035] q _ demand class which indicates relative priority. The 

higher the numerical value of q the lower the relative pri- 
ority. 

[0036] ^_ resource capacity which could be a machine, labor 
hour, or other constraint 

[0037] „ represents a consumer location which refers to an inter- 
nal plant, external demand center, or toa generic indicator 
meaning any plant/or demand center 

[0038] Definition of Objective Function Coefficients 

[0039] PRC - cost of releasing one piece of part m during pe- 

jmae 

riod J at plant a using process e 
[0040] SUBC substitution cost per piece of part number n 

jmna - 

which is being substituted by part number m during pe- 
riods at plant a 

[0041] jQ _ transportation cost per piece of part number m 

jmav 

leaving plant a during period j which are destined for 
plant V 

[0042] iNvc - inventory cost of holding one piece of part 

jma 



number m at the end of period j at a particular plant a 
[0043] DI\/IAXC - cost per piece of exceeding tlie maximum 

jzau 

amount of sliipments of group z parts from plant a to 
consuming location(s) m during period j 
[0044] DMINC - cost per piece of falling short of the minimum 

jzau 

amount of shipments specified for group z parts from 
plant a to consuming location(s) u during period j 
[0045] BOC - backorder cost of one piece of part m at the end 

jmkq 

of period jfor class q demand at customer location k 

[0046] Definition of Constants 

[0047] DEMAND - demand requested during time period j for 

jmkq 

part number m at customer location k for demand class q 
[0048] RECEIPT - quantity of projected wip and purchase order 

jma 

receipts for part number m expected to be received at 
plant during time periodj 
[0049] CAPACITY - Capacity of resource w available at plant a 

jaw 

during period jto support production starts 
[0050] CAPREQ - Capacity of resource w required for part 

Jmaew 

number m at plant a for process e during period jQ 
[0051] TYPER - quantity of component m needed per part 

Jmaen 

number n during period j at plant a using process e 
[0052] YIELD - output of part number m per piece released or 

jmae 

Started at plant a during time period j using process e 



[0053] SUBQTY - quantity of part number m required to sub- 

jmna 

stitute for one piece of part number n at plant a during 
time period j 

[0054] MAXPCT - maximum percentage of total shipments of 

jzau 

group z (collection of parts) leaving supplier a during pe- 
riod jto support consumption at consuming location(s) u 
[0055] MINPCT - minimum percentage of total shipments of 

jzau 

group z (collection of parts) leaving supplier a during pe- 
riod jto support consumption at consuming location(s) u 
[0056] CT - Cycle time. The number of periods between the 

jmae 

release and completion of part m jobs for releases made 
using process e at plant a during time period J 
[0057] jj _ transport time for part number m from plant a to 

mav 

plant V 

[0058] Definition of LP Decision Variables 

[0059] I _ Inventory at the end of period j for part number m at 

jma 

a particular plant a 
[0060] p _ Production starts of part m during period j at plant 

jmae 

a using process e 
[0061] |_ _ Quantity of part number n which is being substi- 

jmna 

tuted by part number m during periodj at plant a 
[0062] J _ Internal shipments of part number m leaving plant 

Jmav 

a during period j which are destined for plant v 



[0063] F - Shipments of part number m leaving plant a dur- 

jmakq 

ing period J and satisfying class ^demand at external cus- 
tomer k 

[0064] B - Back orders of part m at the end of period j for 

jmkq 

class q demand at customer location k 
[0065] H - Total shipments of group z (z is a "collection" of 

jzu 

parts) leaving suppliers during period j to support con- 
sumption at consuming location(s) u. 
[0066] s _ Amount by which total shipments of parts in z from 

jzau 

plant a to consuming location(s) u during period j exceeds 
the maximum amount specified as desired in the sourcing 
rules 

[0067] Q - Amount by which total shipments of group z parts 
from plant a to consuming location(s) u during period j 
falls short of the minimum amount specified as desired in 
the sourcing rules 

[0068] IP Equations or Formulation 

[0069] The following minimizes the objective function subject to 
the constraints shown below. 

[0070] Objective Function: 

[0071] Minimize: 



J m k q 

[0072] Subject to: 

[0073] Sourcing Constraints: 

E + Z ^»«», ) - -S";^ ^ MAXPCTj^ Hj„ 
I " 

m q 
a 

[0074] Capacity Constraints: 
[0075] Backorder Constraints: 

^jn^ = + DEMAND -"^F^^ 

a 

[0076] Material Balance Constraints: 

Ii^ = I(j-iym, + RECEIPTj^+ ^ ^ YIELDxma.* Pm^ + ^ + 
JL 2 T^-T SUBQTrj^*U«.-^ 2] 

isacompama 



[0077] 



Non-Negativity Constraints: 



^, J... ^ 0 , where X is a generic decision variable and /, / etc. represent generic subscripts 

[0078] Other approaches use a heuristic approach to schedule in 
which WIP is sequenced according to its pegged demand 
priority at the various worl<-centers throughout the supply 
chain. While this can be effective from the point of view of 
rigorously considering demand priorities, these heuristics 
typically fail to find globally optimal solutions. From the 
point of view of global optimization of resource allocation, 
linear programming (LP) is a superior approach. However, 
the LP approach relies on cost penalty weightings to "en- 
courage" demand priority compliance and in practice de- 
mand priorities may be violated (e.g. a high yield low pri- 
ority item may be favored over a low yield high priority 
item). 

[0079] The Leachman method (Leachman, R. C, R. F. Benson, C. 
Liu, D.J. Raar. 1996. IMPReSS: An automated production- 
planning and delivery-quotation system at Harris Corpo- 
ration-semiconductor sector. Interfaces 26(1) 637) incor- 
porates the consideration of "preemptive" demand priori- 
ties into an LP modeling framework in which each LP cor- 
responds to a single demand priority. However, this 
method can have severe disadvantages in terms of com- 
putational efficiency, and does not consider hierarchical 



demand prioritizationThe inventive approacli, on tlie otiier 
liand, is a metliod based on aggregating demand priori- 
ties into groups, in wliicli eacli group corresponds to an 
LP run. (Figure 3 illustrates the grouping of priorities into 
sets.) Thus, the present invention is a hybrid method 
which handles hierarchical demand priorities via demand 
class aggregation and subsequent demand tier disaggre- 
gation. The invention generates backorder cost penalties 
associated with demand priorities within a group to model 
the consideration of multiple demand classes within a de- 
mand class group in the LP model. Further, the inventive 
method iteratively modifies and solves the respective LPs 
in memory and leverages the results from one solution to 
the next. 

[0080] The inventive method results in substantial computational 
savings over conventional systems while at the same time 
respecting demand classes throughout the supply chain. 
The invention solves relaxed (e.g., modified) versions of 
the LP model at each iteration that allow for flexibility in 
realigning resources (e.g. work-center capacity) to ac- 
commodate lower priority groups without sacrificing 
higher priority groups. Further, the invention can be ap- 
plied to hierarchical demand priorities by selectively dis- 



aggregating demand tiers associated with demand 
classes. 

[0081] Ideally, the method used for modeling backorder costs 
would guarantee demand compliance for any number of 
demand priorities in a single LP run. In theory this could 
be done by setting the backorder penalty for a particular 
demand priority arbitrarily higher than the penalty for the 
next lower demand priority. However, in practice due to 
finite numerical accuracy of floating point operations 
there are bounds on the allowable range of objective 
function penalties (experimental evidence implies this is 
about .01 to 1,000,000 when operating with double pre- 
cision using the 64 bit C/C++ compilation mode), and the 
delta between priorities which is recognizable. Therefore, 
with the invention backorder costs penalties are deter- 
mined independently for each set (group) of priorities. 
Within each group of demand priorities, the backorder 
penalties comprise a full spectrum range of objective 
function penalties. Thus, each successive linear program- 
ming model allocates a full range of backorder costs 
within the priority group to which the resources are cur- 
rently being allocated to avoid the foregoing problem. 

[0082] The numerical bounds on backorder cost penalties are il- 



lustrated in the example of Figure 2 in whicli tliere are 
lower and upper bounds on the cost penalties based on a 
recursive set of mathematical equations described below. 
Figure 2 further illustrates three possibilities for calibrat- 
ing backorder costs from a continuous range of possibili- 
ties. In Figure 2, the examples illustrate costs which are 
(a) decreasing at an increasing rate with respect to de- 
mand class (b) linearly decreasing in demand class and (c) 
decreasing at a decreasing rate with respect to demand 
class. The inventive modeling backorder costs is based on 
the following recursive equations 
[0083] Delta_{j+l} = C*Delta_{j} 

[0084] where Delta is chosen, such that within each group of pri- 
ority demands, the most and least important demand pri- 
ority correspond to the maximum and minimum allowable 
backorder costs respectively. If there are N demand prior- 
ities in the group and we start with the initial condition 
that backorder cost BC_{l}=max and model the costs as 
follows 

[0085] BC_{1} = MAX_COST 

[0086] BC_{2}=Delta*MAX_C0ST 
[0087] BC_{3}=C*DeltaA{2}*MAX_COST 



[0088] BC_{i} = CA{i-2}*DeltaA{i-l}*MAX_C0ST 

[0089] Since we require tliat BC_{n}=l\/IIN_COST it follows that 

[0090] Delta = 

(MIN_COST/(cA{(n-2)*(n-l)/2}*MAX_COST))A{l/(n-l)} 

[0091] The parameter C is a shape parameter (defined by the 

user) that influences the relative difference between back- 
order costs from one demand priority to the next. Given 
the assumption that backorder costs should be decreasing 
in the importance of the demand priority the values C can 
take are in the following range 

[0092] Lower Bound (LB) = (MIN_COST/MAX_COST)a 
)A{2.0/((n-2)*(n-l))} 

[0093] Upper Bound (UB) = 1/LB 

[0094] where LB and UB follow from the requirement that Delta < 
1 and C*Delta < 1 respectively. Depending on the choice 
for C in the range [LB, UB] different relative differences 
between backorder costs are obtained. These differences 
are illustrated in Figure 2. In the present implementation 
of the invention typically C=LB + 2/3*(UB-LB) is chosen to 
maximize the relative difference between backorder costs 
from one demand priority to the next with some prefer- 
ence (determined by the choice of C) given to the differ- 



ence between more important demand priorities. As C ap- 
proaches UB the preference for relative differences in de- 
mand priority is shifted from less important to more im- 
portant demand classes. 

[0095] The following is a pseudo code description and corre- 
sponding flowchart shown in Figure 5 illustrate an imple- 
mentation of the above described idea. With the invention 
any commercial linear program solver could be employed. 

[0096] Definitions: 

[0097] N(i)=number of demand classes in group i 
[0098] M = number of demand tiers per demand class 
[0099] n = number of demand priority groups 

[0100] F_{mJ,a,k,q}= LP decision variable representing quantity 
of shipments of part m in period j from plant, a, to cus- 
tomer location, v, with priority q 

[0101] B_{m,j,k,q} = LP decision variable representing the quan- 
tity of part m to backorder in period j for customer k of 
priority q 

[0102] D_{m,j,k,q} = demand for part m in period j from cus- 
tomer k of priority q 

[0103] BO = backorder constraints of the form B_{m,j,k,q} = 
B_{m,j-l,k,q} + D_{m,j,k,q} Sum(over a) F_{m,j,a,k,q} 



[0104] BC_{mJ,k,q} = LP cost penalty for backorder variable 
B_{mJ,k,q} 

[0105] MIN.COST = precision based lower bound on BC penalty 
[0106] |\/|AX_C0ST = precision based upper bound on BC penalty 
[0107] UB = upper bound 
[0108] LB = lower bound 

[0109] Step 1: User defines a set of "demand class groups"from 
the set of demand classes. For example, 10 demand 
classes could be split into 2 groups [1,...3] and [4,.., 10]. 
The groups are indexed from i = 1, ...n and the number of 
distinct demand classes in group i is N(i). Set i = 1. 

[0110] Step 2: Fix all F_{mJ,a,k,q} variables where demand class 
q is in a numerically higher (less important) priority group 
than group i and unfix all variables such that q is in a nu- 
merically lower (more important) priority group than 1 or 
equal priority as group 1. Fixing and unfixing is done by 
setting variable_upper_bound = 0.0 and vari- 
able_upper_bound = INFINITY respectively. 

[0111] Step 3: Recalibrate cost penalties for variables, B_{m,j,k,q}, 
with q in group i. This method is based on above discus- 
sion about maximizing the relative difference between 
each priority in group i subject to user defined cost shape 



parameters. 
[0112] //Bounds on shape parameters 

[0113] LB = (MIN_COST/MAX_COST)A{2.0/((N(i)-2)*(N(i)-l))} 
[0114] UB = 1.0/lower_bound 

[0115] //Convert user defined parameter, C, from [0,1] range to 
[LB,UB] 

[0116] Factor = C*LB + (1.0- C)*UB 

[0117] //Compute relative cost difference factor 

[0118] Delta=MIN_COST/(FactorA{(N(i)-2)*(N(j)-l)/2}*MAX_COST) 

)A{1.0/(N(i)-l)}) 
[0119] //Compute penalties from deltas 

[0120] Cost[0] = MAX_COST 

[0121] for(j=l;j<NG);j++) CostGl = factorA{j-i}*Delta*Cost|j-l] 

[0122] //Set BC's to computed cost penalties 

[0123] For(all BC_{m,j,k,q} such that q an element of group i){ 

[0124] DP = relative group i demand priority for current 

BC_{m,j,k,q} 
[0125] BC_{m,j,k,q} = Cost[DP] 

[0126] } 



[0127] Step 4: Solve modified group i LP: 

[0128] if(j= = i) solve LP using normal OSL simplex LP method 

[0^29] else warm start LP solution using previous solution as an 
advanced basis 

[0130] Step 5: Add a new constraint set to the current LP model 
enforcing that B_{mJ,k,q} variables in group i are lower 
bounded based on the current LP solution from Step 4. 
These constraints have the form: 

[0131] For a particular period, j, and demand class q: Sum(parts 
m, customer locations k) B_{mJ,k,q} <= Sum(parts m, 
customer locations k) (B_{mJ,k,q} from Step 4) 

[0132] Step 6: Reset the LP basis with current variables fixed us- 
ing the dual simplex method. Increment group to 1=1+1. If 
(i==n+l) go to step 7, otherwise go to step 2. 

[0133] In the remaining description of the method we have two 
alternative embodiments which we represent as Step 7A 
and Step 7b: 

[0134] Step 7a: For each demand class disaggregate BO con- 
straints and associated variables into separate demand 
tiers. For example, the BO constraint 

[0135] B_{mJ,k,q} = B_{m,j-l,k,q} + D_{m,j,k,q} Sum(over a) 
F_{m,j,a,k,q} 



[0136] would become a set of constraints indexed by tier, t, as 
follows 

[0137] B_{mJ,k,q,t} = B_{mJ-l,k,q,t} + D_{mJ,k,q,t} Sum(over a) 
F-{iTiJ,a,k,q,t} , for all t=l,...,M 

[0138] Thus, the constraints added in Step 5 are disaggregated. 
The previous (aggregated) solution can be used to provide 
a feasible starting solution for reoptimizing the LP with 
the additional disaggregated constraints. For example, if 
there were two demand tiers then F_{m,j,a,k,ql} = Dl* 
F_{m,j,a,k,q} and F_{m,j,a,k,q2} = D2* F_{m,j,a,k,q}, where 
Dl and D2 represent the fraction of demand in period j 
that is in tier 1 and tier 2 respectively. Once the revised 
solution to the disaggregated LP is determined, the F and 
B variables for the current demand class are fixed to their 
current values, the next demand class is disaggregated, 
solved, and so on. This is carried our until all demand 
classes have been considered. 

[0139] Step 7b: For each variable, B_{m,j,k,q}, the invention recal- 
culates penalties for all of the demand classes that re- 
flects priority differences within a class. For example 
rather than add a new demand tier index, t, to variables, 
B-{m,J,k,q}, instead the associated penalty would be a 
weighted average of the different demand tiers for cus- 



tomers, k, within class, q. A final run with these modified 
cost penalties results in a solution which:a) respects de- 
mand class priorities pre-emptively based on the con- 
straints added in Step 5b) given the myriad of solutions 
respecting (a), this final solution will be one which tends 
to respect the differences in demand tierThe alternative 
Steps 7a and 7b allow for the ability to tradeoff different 
criteria. For instance, Step 7a is more accurate in the way 
that it applies demand priorities than Step 7b. On the 
other hand. Step 7b is likely to be computationally more 
efficient since it does not involve the disaggregation of 
variables and constraints that is done in Step 7a. Figure 5 
illustrates the flow of the above steps in the method. 
[0140] Figure 6 is a more general flowchart showing the opera- 
tion of the invention. In item 600, the invention aggre- 
gates the demand priorities into different priority groups. 
In item 602, the invention allocates the resources to the 
highest priority group of demand priorities using a first 
linear programming model. Next, in decision block 604, 
the invention determines if there are additional groups of 
priority demands to process. If so, the invention modifies 
the linear program to create a second linear programming 
model as shown in item 606. Then processing returns to 



item 602 where the invention allocates remaining re- 
sources to the next highest priority group of demand pri- 
orities using a second linear programming model. The 
second linear programming model uses results from the 
first linear programming model. The invention continues 
this process by iteratively repeating the process of allo- 
cating remaining resources to the remaining groups of 
demand priorities in order of priority. After all groups are 
processed (or there are no resources remaining) this por- 
tion of the process is complete 608. 
[0141] When repeating the process of allocating remaining re- 
sources, the method uses a different linear programming 
model for each iteration and each different linear pro- 
gramming model in item 606 uses results of the previous 
linear programming model. In each of the linear program- 
ming models created in item 606, the invention fixes vari- 
ables associated with all demand priorities that are in 
lower priority (numerically higher) groups than the priority 
group to which the resources are currently being allo- 
cated. Also, each linear programming model in item 606 
allocates the full range of backorder costs within the sin- 
gle priority group to which the resources are currently be- 
ing allocated. Therefore, each iteration has the benefit of 



a full range of backorder costs. The invention can also di- 
vide the priority groups Into different sub-prlorlty tiers. 
These sub-prlorlty tiers can be processed simultaneously 
or separately. 

[0142] The present Invention has been implemented on an IBM 

P690 server using the AIX operating system. The steps for 
implementing the present invention are preferably pro- 
grammed in C/C+ + . It should be understood by those of 
ordinary skill in the art, however, that the represent in- 
vention is not limited to the above implementation and is 
independent of the computer/system architecture. Ac- 
cordingly, the present invention may equally be imple- 
mented on other computing platforms, programming lan- 
guages and operating systems, and also may be hard- 
wired Into a circuit or other computational component. 

[0143] The present Invention Is a hybrid method which handles 
hierarchical demand priorities via demand class aggrega- 
tion and subsequent demand tier dIsaggregatlon.The in- 
ventive approach Is based on aggregating demand priori- 
ties Into groups, In which each group corresponds to an 
LP run. The Invention generates backorder cost penalties 
associated with demand priorities within a group to model 
the consideration of multiple demand classes within a de- 



mand class group in the LP model. Further, the invention 
method iteratively modifies and solves the respective LPs 
in memory and leverages the results from one solution to 
the next. 

[0144] jhe inventive method results in substantial computational 
savings over conventional systems while at the same time 
respecting demand classes throughout the supply chain. 
The invention solves relaxed (e.g., modified) versions of 
the LP model at each iteration that allow for flexibility in 
realigning resources (e.g. work-center capacity) to ac- 
commodate lower priority groups without sacrificing 
higher priority groups. Further, the invention can be ap- 
plied to hierarchical demand priorities by selectively dis- 
aggregating demand tiers associated with demand 
classes. 

[0145] While the invention has been described in terms of the 
preferred embodiments, those skilled in the art will rec- 
ognize that the invention can be practiced with modifica- 
tion within the spirit and scope of the appended claims. 



